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Summary 


The  objective  of  this  study  is  determination  of  the  accuracy 
with  which  the  station-source  backazimuth  can  be  determined  from 
single  station,  three-component  data  in  the  case  of  industrial 
explosions  at  near  regional  distance.  The  single  station  backazimuth, 
derived  from  the  polarization  of  the  P  wave,  was  a  parameter 
reported  by  U.S.  stations  during  the  GSETT-2  experiment.  One  of 
these  stations  (BLA),  located  at  Blacksburg,  Virginia,  is  situated  to 
record  numerous  large  mining  explosions  in  the  distance  range  100 
to  300  km.  The  explosions  can  be  reliably  located  using  arrival  time 
data  from  a  regional  network. 

A  moving  window  polarization  analysis  of  the  explosion  P  wave 
arrivals  was  performed.  In  a  fashion  similar  to  the  procedures 
followed  during  the  GSETT-2  experiment,  maximum  rectilinearity  of 
motion  was  used  as  a  criterion  for  judging  the  "best"  backazimuth 
estimate  from  the  three-component  P  wave  data.  The  resulting  best 
estimate  backazimuth  is  compared  with  that  derived  from  the 
independently  located  mines. 

The  results  of  the  comparisons,  performed  to  date  on  37 
signals,  show  that  three-component  signal/noise  ratios  less  than  2.0 
can  result  in  very  unreliable  backazimuths.  However,  for  35  signals 
with  S/N  ratios  exceeding  2.0,  the  mean  backazimuth  error  is  6 
degrees  with  standard  deviation  21  degrees.  The  error  in  the 
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determinations  does  not  appear  to  be  strongly  correlated  with 
signal/noise  ratios  in  the  range  2  to  10. 


Characteristics  of  the  signals  which  contribute  to  the  errors  in 
the  single-station  backazimuth  estimates  include: 

1.  Very  emergent  initial  motions  from  the  delay-fired 
explosions. 

2.  Off-azimuth  P  wave  arrivals  and/or  converted  phases 
arriving  very  early  in  the  P  wave  coda. 

3.  Steep  apparent  angles  of  incidence,  averaging  22  degrees 
from  vertical,  which  tend  to  reduce  the  signal/noise  ratios 
on  the  horizontal  components. 

Optimum  data  segments  for  reliable  source  backazimuth 
estimates  are  restricted  to  very  short  time  intervals  (1  second  or 
less)  beginning  with  the  initial  P  wave  motion. 


1. 


Introduction 


As  part  of  the  GSETT-2  experiment,  a  high-quality  six-channel 
(three  component  short-period  and  broadband)  seismic  system  was 
installed  at  Blacksburg,  Virginia,  in  January,  1990.  The  Blacksburg 
location  is  well  situated  to  record  industrial  explosions  in  the 

distance  range  100  to  300  km  arising  from  surface  coal  mining 
operations  in  Ohio,  West  Virginia  and  Kentucky. 

This  report  examines  the  accuracy  with  which  the  source- 
station  backazimuth  can  be  determined  at  station  BLA  from 
polarization  analysis  of  the  mining  explosion  generated  P-wave 

signals.  As  a  basis  for  the  accuracy  assessment,  the  backazimuths 

derived  from  the  three  component  analysis  at  BLA  are  compared 
with  the  locations  of  the  mines  estimated  using  arrival  time  data 
recorded  at  stations  of  the  Virginia  Tech  and  Tennessee  Valley 

Authority  regional  seismic  networks. 

2.  The  Blacksburg  Three-Component  Station 

The  Blacksburg,  Virginia,  GSETT-2  station  was  installed  in 
January,  1990  in  the  WWSSN  vault  on  the  campus  of  Virginia  Tech. 
This  vault  is  founded  on  Cambrian  dolomite  of  the  Appalachian 
Valley  and  Ridge  geologic  province.  The  rtatinn  is  <*qiiipped  with 
Teledyne  Geotech  GS-13  short-period  and  BB-13  broadband  sensors. 
Digitization,  multiplexing  and  calibration  are  performed  by  a 
Teledyne-Geotech  RDAS-200  unit,  and  time  synchronization  is 
obtained  via  a  Kinemetrics,  Inc.  Omega  rcceiver/clock.  The  short- 
period  and  broadband  sample  rates  are  40  and  10  samples/sec. 
respectively.  The  data  are  demultiplexed,  event-detected  and 
archived  by  a  Science  Horizons,  Inc.  NOMAD  workstation. 

2.1  Calibration 

Relative  differences  in  system  response  on  the  two  horizontal 
components  result  in  systematic  errors  in  backazimuth  estimates. 


Hence,  the  system  responses  of  the  short-period  channels  were 
determined. 

The  RDAS-200  unit  provides  for  the  capability  of  either  pulse- 
type  transient  or  sinusoidal  steady  state  calibration  test  signals  to  be 
applied  to  the  sensor  calibration  circuits.  The  transient  calibration 
pulses  were  used  to  determine  the  short-period  system  response 
using  the  approach  described  by  Chapman  et  al.  (1988).  Essentially, 
the  method  fits  an  ideal  damped  oscillator  amplitude  response  to  the 
Fourier  amplitude  spectrum  of  a  recorded  system  calibration  pulse. 
A  nonlinear  least  squares  algorithm  is  used  to  estimate  the  system 
period,  damping  and  gain. 

The  amplitude  responses  of  the  horizontal  channels  were  found 
to  differ  by  less  than  2  percent  in  the  1  to  10  Hz  frequency  band. 
The  vertical  component  response  differs  from  that  of  the  two 
horizontals  by  approximately  9  percent.  These  differences  were 
removed  by  mathematical  correction  of  the  recorded  data. 

3.  The  Regional  Network 

The  network  stations  used  to  determine  the  locations  of  mining 
explosions  are  shown  in  Figure  1.  The  six  stations  in  western 
Virginia  and  in  West  Virginia  are  operated  by  Virginia  Tech,  whereas 
the  stations  in  Tennessee  and  Kentucky  are  part  of  a  network 

operated  by  the  Tennessee  Valley  Authority  (TVA).  The  Virginia 
Tech  stations  are  recorded  digitally  in  an  event  triggered  mode.  The 

TVA  stations  are  recorded  on  16  mm  Develocorder  photographic  film. 

4.  Location  and  Characteristics  of  the  Mining  Explosions 

Figure  1  shows  the  locations  of  37  mining  explosions  processed 
to  date.  Examination  of  the  signals  using  trace  overlays  recorded  by 
the  Virginia  Tech  network  indicates  that  the  data  set  consists  of 

multiple  explosions  at  nine  mine  sites  plus  single  explosions  at  nine 
additional  sites.  Confirmation  of  the  geographic  locations  of  the 


explosions  through  contact  with  the  mine  operators  has  been 
achieved  tor  three  mines.  These  three  confirmed  locations  are  the 
Central  Ohio  Mine  ia  southeastern  Ohio,  the  Ruffner  .Mine  in  West 
Virginia  and  the  -.-tarfire  Mine  in  eastern  Kentucky,  shown  in  Figure 

n 

The  epicenters  of  the  explosions  were  located  using  the 
c,.(mbined  Virginia  Tech  and  TVA  network  arrival  time  data  (both  P 
and  S-Lg  arrival  times).  The  location  program  HYPOELLIPSE  (Lahr, 
1980)  veas  used  with  a  three  layer  crustal  velocity  model  (Bollinger 
et  al.,  1980).  The  network  locations  and  location  quality  estimates 
are  given  in  Table  1.  In  cases  where  the  event  epicenter  was  not 
known  from  independent  information  and  where  multiple  events 
were  recorded  from  the  same  location,  the  network  location  with 
smallest  statistical  uncertainty  was  adopted  as  the  epicenter  of  the 
group  of  explosions. 


rc  I  rjC*? 


Figure  1.  Network  stations  (including  three-component  station  BLAi 
are  shown  by  circles.  Crosses  indicate  mine  locations 
yielding  multiple  explosion  signals.  Mines  for  Ahich  onlv 
a  single  explosion  was  recorded  are  shown  by  tii.mgles. 


Table  1 


EVENT 

LAT 

LONG 

ERH2 

DIST 

BAZ 

BAZ 

NO. 

Mine^ 

fDes.  N) 

mep.  W)  rkm^ 

Ckm) 

(Deg.) 

fEst) 

RECT 

s.x 

20 

A 

39.80 

81.69 

308 

339 

334 

0.88 

1 . 1 

2  1 

A 

39.80 

81.69 

- 

308 

339 

368 

0.91 

1.4 

3  1 

A 

39.80 

81.69 

- 

308 

339 

20 

0.93 

1.7 

35 

A 

39.80 

81.69 

- 

308 

339 

239 

0.87 

1 .7 

36 

A 

39.80 

81.69 

- 

308 

339 

33  1 

0.89 

1  0 

B 

37.43 

83.12 

- 

241 

277 

275 

0.99 

3.U 

1  3 

B 

37.43 

83.12 

- 

241 

277 

250 

0.98 

7  ^ 

1  6 

B 

37.43 

83.12 

- 

241 

277 

270 

0.96 

7  ^ 

28 

B 

37.43 

83.12 

- 

241 

277 

289 

0.97 

4.1 

44 

B 

37.43 

83.12 

- 

241 

277 

243 

0.97 

3.S 

1  5 

C 

37.80 

81.81 

- 

139 

298 

248 

0.94 

6, 1 

47 

C 

37.80 

81.81 

- 

139 

298 

27  1 

0.98 

3.S 

1  9 

- 

38.53 

80.39 

1  5 

146 

1±4,2 

357 

0.88 

1 .8 

22 

- 

38.53 

80.39 

1  1 

146 

1±4.2 

9 

0.85 

3.2 

1  2 

- 

38.13 

82.40 

7 

202 

301±1.5 

299 

0.98 

6.9 

3  9 

- 

38.13 

82.40 

7 

202 

301±1.5 

286 

0.98 

5.0 

43 

- 

38.13 

82.40 

2  1 

202 

301±1.5 

313 

0.94 

3 . 1 

I  8 

- 

37.99 

82.54 

5 

206 

295±0.7 

280 

0.97 

'I  ^ 

3  2 

- 

37.99 

82.54 

3 

206 

295±0.7 

265 

0.98 

9 

24 

- 

37.58 

82.18 

8 

1  6  1 

285±2.8 

144 

0.93 

1  .4 

3  8 

- 

37.58 

82.18 

9 

1  6  1 

285±2.8 

283 

0.93 

1.3 

4  1 

37.58 

82.18 

1  0 

161 

285±2.8 

3  1  4 

0.96 

5 .4 
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Table  1  (continued) 


EVENT 

LAT 

LONG 

EREP 

DIST 

BAZ 

BAZ 

NO 

M  i  n  e  ^ 

(Deo.  N)  (Deo.  W)  (kmi 

(km) 

(Dec. ) 

(E,st) 

RFCT 

S  N 

1  7 

38.03 

81.69 

9 

1  44 

309±i .6 

307 

0.96 

1  1 

42 

- 

38.03 

81.69 

9 

1  44 

309±i.6 

3  I  1 

0.91 

2,5 

34 

- 

38.03 

81.69 

5 

1  44 

309±1 .6 

320 

0.90 

2,5 

2  7 

- 

37.96 

81.88 

1  4 

154 

303-2.1 

287 

0.94 

3.0 

:>  j 

- 

37.96 

81.88 

/ 

154 

303±2.I 

308 

0.97 

8.7 

8 

- 

38.34 

82.78 

1  7 

243 

302±3.6 

308 

0.92 

2.0 

9 

- 

37.52 

82.64 

99 

200 

280±2.3 

279 

0.96 

4.8 

1  1 

- 

38.04 

82.26 

1  1 

1  87 

300-2.2 

273 

0.95 

3,8 

1  4 

- 

37.94 

82.65 

6 

2  1  2 

293±1.1 

264 

0.96 

6.9 

25 

- 

37.19 

81.51 

1  2 

97 

269±4.5 

274 

0.94 

•i 

26 

- 

37.85 

82.54 

7 

200 

292±1 .7 

284 

0.90 

*« 

29 

- 

38.28 

82.36 

4 

208 

305±1.0 

279 

0.97 

".9 

37 

- 

38.28 

81.21 

6 

137 

330±2.2 

336 

0.S9 

^  1 

40 

- 

38.04 

82.16 

1  9 

179 

304±3.7 

■*»  5  ? 

0.9~ 

-1 

•i 

lA  = 

Central  Ohio  Coal  Company; 

B  = 

"Starfire 

Mine, 

eastern  Kentucky: 

C  = 

"Ruffner’ 

Mine, 

Arch  Minerals,  Inc. 

,,  Yolyn,  West  Vir, 

ginia; 

-  =  Mine  location  derived  from  network  arrival  time  data:  not  confirmed 
by  independent  information 

^ERH  =  Horizontal  error  measure  from  location  program  HYPOELLIPSE. 

^Errors  estimated  on  the  basis  of  the  angles  subtended  by  the  HYPOELLIPSE 
68%  confidence  error  ellipses. 
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The  location  capabilities  of  the  network  were  examined  by 
comparing  the  network  derived  epicenter  locations  of  the  Starfire. 
Ruffner  and  Central  Ohio  Mines  with  their  known  locations.  Figure  2 
shows  the  HYPOELLIPSE  68%  chi-square  confidence  regions  f<  the 
horizontal  projections  of  the  hypocenter  locations.  Also  shown  in  the 
figure  are  the  actual  locations  of  the  three  mines.  The  network 
locations  for  the  explosions  at  the  Starfire  Mine  in  Kentucky  and  the 
explosions  at  the  Ruffner  Mine  in  West  Virginia  are  in  good 
agreement  with  the  actual  locations.  The  network  location 'capability 
for  the  Ohio  explosions  is  degraded  because  of  very  low  signal-to- 
noise  ratios  for  that  more  distant  source. 

Figure  3  shows  the  location  error  ellipses  for  the  remaining 
events  in  the  data  set.  With  the  exception  of  Event  9  in  Kentucky, 
which  was  well  recorded  only  by  the  Virginia  Tech  network,  the 
error  ellipse  semi-major  axes  are  all  less  than  25  km.  The 
uncertainty  in  the  network  derived  backazimuths  must  be  taken  into 
consideration  in  order  to  assess  the  quality  of  the  single-station 
estimates.  The  horizontal  angles  subtended  by  the  HYPOELLIPSE 
error  ellipses  as  viewed  from  station  BLA  were  calculated  and  are 
listed  in  Table  1.  The  maximum  uncertainty,  as  defined  by  the 
angular  extent  of  the  error  ellipse,  is  ±4.6  degrees,  for  Event  #25. 
This  error  estimate  does  not  account  for  any  systematic  bias  due  to 
inaccuracy  of  the  velocity  model.  However,  this  latter  source  of  error 
appears  to  be  small,  for  as  shown  in  Figure  2,  the  calculated  locations 
of  the  Starfire  and  Ruffner  Mines  are  near  the  true  mine  locations. 

Information  obtained  from  the  mine  operators  at  the  three 
mines  indicates  that  the  larger  explosions  usually  are  the  result  of 
"cast  blasting,"  a  technique  which  typically  utilizes  delay-fired, 
multi-row  charge  designs,  similar  to  the  explosions  studied  by 
Chapman  et  al.  (1992). 
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Confirmed  mine  locations  are  shown  bv  the  circles. 

¥ 

Location  of  station  BLA  shown  by  the  triangle.  Ellipses 
indicate  the  horizontal  projections  of  the  68%  chi-square 
confidence  ellipsoids  for  blast  epicenters  from  location 
program  HYPOELLIPSE. 
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Confirmed  mine  locations  are  shown  by  the  circles  with 
associated  letters  as  listed  in  Table  1.  The  location  of  the 
three-component  station  is  shown  by  the  triangle.  The 
ellipses  indicate  the  horizontal  projections  of  the 
HYPOELLIPSE  68%  chi-square  confidence  ellipsoids  for 
mine  locations  inferred  from  blast  epicenters. 
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Method  of  Polarization  Analysis 


The  backazimuth  estimates  from  the  three-component  BLA 
station  were  derived  using  the  method  developed  by  Jurkevics 
(1988).  The  processing  is  carried  out  in  the  time  domain  using  a 
series  of  overlapping  time  windows.  For  each  time  window,  a  3  x  3 
data  covariance  matrix  is  formed  from  the  auto-  and  cross -variances 
of  the  three  components  of  motion.  The  principle  axes  of  the 
polarization  ellipsoid,  (Xj  uj,  i  =  1,  2,  3)  for  each  window  are 
determined  by  solving  the  eigenproblem  for  eigenvalues  (A.i,X2and 
X2))  and  eigenvectors  (ui,  U2,  and  U3).  For  P  wave  motion,  the  azimuth 
of  propagation  is  estimated  from  the  horizontal  orientation  of  the 
eigenvector  corresponding  to  the  largest  eigenvalue,  ui: 

Pazimuth  =  tan'l  (U2i/U3i),  (1) 

where  uji,  j  =1,2,3  are  the  direction  cosines  of  the  eigenvector 
associated  with  the  largest  eigenvalue.  A  measure  of  the  degree  of 
rectilinearity  of  the  particle  motion  is  given  by 

R=l-((X^+i,)/2X,).  (2) 


which  is  1.0  when  X2  and  X.3  are  zero,  for  pure  P  wave  motion.  The 
method  can  be  made  frequency  dependent  by  applying  the  analysis 
to  narrow  bandpassed  time  series  data. 
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6.  Data  Analysis 


The  BLA  short-period  data  were  corrected  for  minor 
differences  in  system  response  of  the  three  component  channels. 
Frequency  domain  assessment  of  the  P  wave  to  pre-P  wave 
signal/noise  levels  was  made  by  calculating  the  Fourier  amplitude 
spectrum  of  the  P  wave  and  pre-P  wave  data  segments  using  10  sec 
time  windows.  Explosions  with  signal/noise  ratios  less  than  2  in  the 
frequency  range  1  to  10  Hz  were  rejected.  The  three  component  data 
were  then  bandpass  filtered  using  a  four-pole  Butterworth  filter  with 
comer  frequencies  1.0  and  10.0  Hz. 

The  duration  of  the  moving  time  window  used  in  the 
polarization  analysis  is  a  critical  parameter  for  determining  accurate 
backazimuth  estimates.  Its  value  will  depend  upon  the  bandwidth  of 
the  signal,  the  sample  rate  and  the  polarization  characteristics  of  the 
signal.  Trial  and  error  testing  showed  that  a  window  duration  of  0.5 
seconds  generally  provided  the  best  backazimuth  determination  with 
our  data  set.  Given  the  data  bandwidth  of  1  to  10  Hz  and  a  sample 
rate  of  40  samples/second,  windows  less  that  0.5  seconds  in  duration 
gave  unstable  results  due  to  the  small  number  of  sample  points, 
whereas  larger  windows  tended  to  reduce  resolution  because  of  the 
non-stationary  polarization  character  of  the  coda  immediately 
following  the  initial  P  wave  arrival. 

The  polarization  analysis  was  performed  in  an  automated 
manner,  with  the  best-estimate  of  the  backazimuth  determined  on 
the  basis  of  maximum  rectilinearity  near  the  time  of  initial  P  wave 
motion. 
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Examples  of  the  various  steps  involved  in  the  processing 
procedure  are  shown  in  Figures  4  through  6.  Shown  in  Figure  4  are 
the  data  from  Explosion  #16  at  the  Starfire  Mine,  rotated  into 
vertical,  radial  and  transverse  components.  Figure  5  shows  the 
Fourier  amplitude  spectra  of  the  10  second  P  wave  and  pre-P  wave 
noise  segments.  Figure  6  shows  a  5  second  duration  sequence  of  the 
data  centered  near  the  time  of  P  wave  onset,  with  continuous 
running  estimates  of  backazimuth,  rectilinearity  and  signal/noise 
ratio.  The  signal/noise  quantity  was  calculated  as  the  ratio  of  the 
average  three-component  amplitude  in  the  5  second  moving  window 
to  the  average  long-term  three-component  amplitude  using  a  10 
second  data  segment  prior  to  P  wave  onset. 
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Figure  4.  Time  series  data  for  Explosion  #16  at  the  Starfire  Mine  in 
Kentucky.  The  data  have  been  rotated  into  vertical  (top), 
radial  (center),  and  transverse  (bottom)  components. 
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25-MAR-92  20;  13:39.00  BLA  sz 
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Fourier  amplitude  spectra  of  the  vertical  component  P  wave 
arrival.  Explosion  #16  (solid  line),  and  the  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  6.  Results  of  polarization  processing  for  Explosion  #16.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial.  (3)  transverse.  {4.) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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7,  Results 


The  backazimuth  errors,  defined  as  the  differences  in  the 
network  derived  estimates  and  the  results  from  the  three- 
component  polarization  measurements,  are  plotted  as  functions  of 
signal/noise  ratio  and  rectiiinearity  in  Figures  7  and  8,  respectively. 
Figure  7  suggests  that  ratios  less  than  2  can,  but  do  not  always,  result 
in  very  unstable  estimates  of  the  polarization  parameters,  as 
evidenced  by  the  two  points  with  backazimuth  errors  in  excess  of  90 
degrees.  These  two  events  are  #35  and  #24  (Table  1).  Note  that 
other  examples  of  explosions  at  both  of  these  mines,  featuring  similar 
signal/noise  ratios,  produced  backazimuth  estimates  much  closer  to 
the  true  values.  An  example  involves  Events  35  and  20,  both  at  the 
Central  Ohio  Mine.  Figures  9  and  10  show  the  three-component 
recordings  for  the  explosions.  Figures  11  and  12  show  the  vertical 
component  P  wave  and  pre-P  wave  noise  spectra  of  the  two  signals, 
and  Figures  13  and  14  summarize  the  three-component  processing 
results  for  both  events.  Note  that  the  moving  window  backazimuth 
estimates  near  the  times  of  P  wave  onset  are  stable  at  values  near 
the  true  value  of  180  degrees  for  the  rotated  components,  in  the  case 
of  Event  20  (Figure  14),  whereas  those  for  Event  35  (Figure  13) 
change  rapidly  at  times  near  the  P  wave  onset.  Clearly,  the  random 
changes  in  the  character  of  the  noise  may  severely  affect  the 
backazimuth  determination,  when  the  signal/noise  ratio  ic  small. 

Large  signal/noise  ratios  do  not  guarantee  accurate 
backazimuth  estimates,  at  least  for  BLA  and  this  mine  explosion  data 
set.  Figure  7  shows  no  evidence  for  the  backazimuth  error  to 
decrease  with  increasing  signal/noise  ratio,  in  the  range  2  to  10. 
Deleting  the  two  observations  with  errors  in  excess  of  100  degrees 
results  in  a  mean  azimuth  error  of  -6  degrees;  the  standard  deviation 
is  21  degrees. 
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Signal/Noise 


Figure  7.  Three-component  station  backazimuth  error  versus 
signal/noise  ratio. 


Rectilinearity 


Figure  8.  Three-component  station  backazimuth  error  versus 
rectilinearity: 
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Figure  9.  Time  series  data  for  Explosions  #35  at  the  Central  Ohio 
Mine.  (Top)  Vertical  component,  (middle)  radial 
component,  (bottom)  transverse  component. 
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Figure  10.  Time  series  data  for  Explosion  #20  at  the  Central  Ohio 
Mine.  (Top)  Vertical  component,  (middle)  radial 
component,  (bottom)  transverse  component. 
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Figure  11.  Fourier  amplitude  spectrum.  Explosion  #35:  Vertical 

component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  12.  Fourier  amplitude  spectrum.  Explosion  #20:  Vertical 

component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  13.  Results  of  polarization  processing  for  Explosion  #35.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4-) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  vvhich  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  14.  Results  of  polarization  processing  for  E.xplosion  #20.  From  the 
top,  the  traces  are;  (I)  Vertical,  (2)  radial.  (3)  transverse.  (4) 
station-source  backazimuth  (deg).  (5)  rectilinearity.  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  8  shows  the  backazimuth  error  plotted  as  a  function  of 

rectilinearity.  Again,  there  is  no  tendency  in  this  data  set  for  the 

errors  to  decrease  with  increasing  values  of  this  parameter  in  the 
range  0.8  to  1.0.  Thus,  it  appears  that  the  P  wave  arrivals  from  some 
of  the  mine  explosions  are  arriving  from  off-azimuth  directions, 
possibly  due  to  refraction  from  laterally  heterogeneous  velocity 
structure. 

As  a  rule,  the  explosion  data  set  exhibits  very  emergent  and 
complex  P  wave  arrivals.  The  emergent  character  of  the  signal  may 
be  due  to  the  extended  nature  of  the  delay-fired  source-time 
function.  For  example,  typical  large  explosions  at  the  Starfire  Mine 
involve  up  to  60  separate  charges,  and  the  duration  of  the  source 
time  function  can  be  in  excess  of  1  second  (Chapman,  et  al.,  1992). 

This  is  also  the  case  for  some  of  the  larger  explosions  at  the  Ruffner 

Mine  (personal  communication,  mine  engineers)  and  is  probably 
typical  of  most  of  the  explosions  in  our  data  set  from  other  mines  as 
well.  Results  from  the  polarization  analysis  indicate  that,  in  general, 
only  the  initial  portion  of  the  P  wave  arrival  exhibits  strongly 
rectilinear  motion,  which  is  polarized  in  the  source-station  azimuth. 
Off-azimuth  arrivals  consisting  of  converted  and  scattered  energy 
appear  at  a  very  early  stage  in  the  coda:  rarely  can  accurate 
estimates  of  the  P  wave  backazimuth  be  obtained  from  portions  of 
the  coda  at  times  greater  than  1.0  second  following  the  initial  onset 
time.  Also,  as  shown  in  Figure  15,  the  P  wave  angles  of  incidence 
with  the  vertical  are  small,  averaging  22  degrees  with  a  standard 
deviation  of  9  degrees.  This  results  in  a  reduced  signal/r.oise  ratio 
on  the  horizontal  components,  which  is  approximately  one-half  of  the 
total  three-component  value  given  in  Table  1.  This  phenomenon, 
which  may  be  station  dependent,  tends  to  further  reduce  the 
accuracy  of  the  backazimuth  determination  for  low  S/N  events. 
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Figure  15.  Apparent  angle  of  incident  P  wave  versus  three- 
component  signal/noise  ratio. 

The  above  complications  make  accurate  single-station 
backazimuth  determination  a  difficult  problem  for  this  delay-fired 
explosion  data  set.  Figures  16  through  27  show  live  examples  with 
signal/noise  ratios  exceeding  4  wherein  the  polarization  analysis 
resulted  in  backazimuth  estimates  correct  to  within  10  degrees  of  the 
true  value.  Note  that  in  most  cases,  the  very  early  portion  of  the 
signal  yielded  the  best  estimate  of  the  backazimuth,  and  that  the 
polarization  characteristics  of  the  immediate  P  wave  coda  tend  to 
fluctuate  strongly. 
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Figure  18.  Results  of  polarization  processing  for  Explosion  #9.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimulh  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  19.  Time  series  data  for  Explosion  #12.  (Top)  Vertical  component, 
(middle)  radial,  (bottom)  transverse. 
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Figure  20.  Fourier  amplitude  spectrum.  Explosion  #12:  Vertical 

component  P  wave  arrival  (solid)  and  pre-P  wave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  21.  Results  of  polarization  processing  for  Explosion  #12.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse.  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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Figure  22.  Time  series  data  for  Explosion  #17,  (Top)  Vertical  component, 
(middle)  radial,  (bottom)  transverse. 
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Figure  24.  Results  of  polarization  processing  for  Explosion  #17.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 


30 


40 


60 


80 


Time  sec 


Figure  25.  Time  series  data  for  Explosion  #33.  (Top)  Vertical  component, 
(middle)  radial,  (bottom)  transverse. 
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Figure  26.  Fourier  amplitude  spectrum.  Explosion  #33:  Vertical 

component  P  wave  arrival  (solid)  and  pre-P  w-ave  noise 
(dashed).  Ten  second  time  windows  were  used. 
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Figure  27.  Results  of  polarization  processing  for  Explosion  #33.  From  the 
top,  the  traces  are:  (1)  Vertical,  (2)  radial,  (3)  transverse,  (4) 
station-source  backazimuth  (deg),  (5)  rectilinearity,  and  (6) 
three-component  signal/noise  ratio.  Note  that  for  rotated 
components,  the  true  backazimuth  angle  is  180  degrees. 

Vertical  lines  indicate  time  at  which  the  rectilinearity  is 
maximum.  The  duration  of  the  moving  window  for  polarization 
analysis  is  0.5  seconds. 
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8.  Conclusions 

We  have  compared  the  single-station  backazimuth  estimates  at 
BLA  with  independently  derived  locations  for  a  set  of  mining 
explosions  and  find  that  for  sources  in  the  northwest  quadrant,  • 

polarization  analysis  of  the  P  wave  arrival  gives  an  essentially 
unbiased  estimate  of  the  source  direction  (mean  error  6  degrees),  but 
that  the  scatter  in  the  results  is  substantial,  with  a  standard 
deviation  for  35  measurements  being  21  degrees. 

Characteristics  of  the  signals  which  contribute  to  the  errors  in 
the  single-station  backazimuth  estimates  include: 

1.  Very  emergent  initial  motions  from  the  delay-fired 
explosions. 

2.  Off-azimuth  P  wave  arrivals  and/or  converted  phases 
arriving  very  early  in  the  P  wave  coda. 

3.  Steep  angles  of  incidence,  averaging  22  degrees  from 
vertical,  which  tend  to  reduce  the  signal/noise  ratios  on 
the  horizontal  components. 

Optimum  data  segments  for  reliable  source  backazimuth 
estimates  are  restricted  to  very  short  time  intervals  (1  second  or 

less)  beginning  with  the  initial  P  wave  motion. 
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